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The results of approximate N 3 LO and detailed NNLO fits to xF 3 data of the CCFR’97 collaboration are pre¬ 
sented. We demonstrate that 1/Q 2 non-perturbative corrections to xF 3 modeled by three independent procedures 
are shadowed by perturbative QCD effects, starting at the NNLO. Special attention is paid to revealing the role 
of the recently calculated NNLO corrections to the anomalous dimensions and N' ! LO corrections to the coefficient 
functions of odd moments of xF 3 with n < 13. The related values of a s (Mz) are extracted. 


It is known that the leading non-perturbative 
power suppressed corrections to DIS structure 
functions (SFs) have the dimension 1/Q 2 . How¬ 
ever, it turned out that phenomenological value 
of non-perturbative effects depend crucially from 
the order of the corresponding perturbative con¬ 
tributions. It is worth to remind that in 1979, 
when the data for DIS neutrino-nucleon scat¬ 
tering was not precise enough, the authors of 
Ref. [ Q] were unable to separate perturbative 
l/ln(Q 2 ) source of scaling violation from the 
l/Q 2 -effects. At present both the precision of 
xF 3 measurements [ |) and the information on 
renormalization-group perturbative QCD evolu¬ 
tion of Mellin moments became more precise. The 
latter ones enter into the Jacobi polynomial for¬ 
mula [ Q] 
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where u>(a,/3) = x a (l — x)P , c^\a,f3) contain 
Euler T functions from a and /3 and the moments 

M™ C {Q 2 ) = M n (Q 2 ) + ( 2 ) 

are taking into account the leading order target 
mass corrections. The kinematic contributions of 
order 1/Q 4 did not affect the results of our previ¬ 
ous less detailed NNLO fits to CCFR’97 data (see 
Refs.[ |], |) and the analysis of Ref.[ |) described 
below. 

In this talk we will concentrate on the results of 
the most recent analysis of the CCFR’97 data for 
xF 3 in the NLO, NNLO and approximate N 3 LO 
levels of perturbative QCD [ || , paying special at¬ 
tention to the possibility of the detection of non- 
perturbative l/Q 2 -contributions to xF 3 . They 
will be modeled by three independent ways. First, 
is the infrared-renormalon (IRR) model of Ref. [ 

0 

TTHH Nmax FI 

= w (a,P) ©1 Q,/3) Z c j"W) M /£?(3) 

n=0 j =0 

where M/ rr = CM n (Q 2 )A' 2 /Q 2 and C = —n + 
4 + 2/(n + 1) + 4 (n + 2) + 4Si(n), calculated in 
Ref. [ fij from the single chain of quark loop inser- 
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tions into the one-gluon contribution to the cor¬ 
responding Born diagram and A 2 is the arbitrary 
fitted parameter. Next, following the NNLO 
Bernstein polynomial fits to CCFR’97 data of 
Ref. [ H (see also Ref. [ 0 ) we consider gradient 
model of twist-4 term, namely 

M«J F 3 (Q 2 )=n^M n (Q 2 ) (4) 

where B 2 is the free parameter. Another possi¬ 
bility is to choose 


HT = h{x ) in the model — independent way .(5) 


Here h{x) is defined by free parameters hi = 
h(xi ), where x, t are the points in experimen¬ 
tal data binning. In our work the following 
renormalization-group equation for the Mellin 
moments of xF, 3 was used: 


M n (Q 2 ) 


M n {Ql) 


2T = eX P 


7 ff(*) 

P(t) 


dt. 


C^{A S {Q 2 )) 


C^(MQD) 


( 6 ) 


where A s = a s /( 47 r) is the MS -scheme coupling 
constant and M n (Q q) is defined in the initial scale 
as M n {Ql) = x n -' 2 A(Ql)x b( - Q °\l+"/{Ql))dx 

At the N 3 LO the expression for C^(A S ) can 
be presented in the following form 

C ( pJ = 1 + C (1) {n)A s + C (2) (n)A 2 + C (3) (n) A 3 (7) 

where the NNLO correction C^ 2 \n) can be ob¬ 
tained for any n from the results of Ref. [ |l0| , 
which were checked with the help of other meth¬ 
ods in Ref. [ [lT[ | . The N 3 LO contributions to Eq. 


were analytically calculated in Ref. [ 12 for 


odd n < 13. The N 3 LO expansion of the anoma¬ 
lous dimension term has the following form 


exp 


f A s(Q 2 ) 7 W( t ) 


/3(f) 


= (A s (Q 2 )) 7 ™ //3 ° x 


dt. 


1 +p(n)A s (Q 2 ) 


+q(n)A s (Q 2 ) 2 + r(n)A s (Q 2 ) 3 


( 8 ) 


where p(n ), q(n) and r(n) are defined through 
coefficients of QCD /3-function and anomalous di¬ 
mension 7 ^ (see Ref. [ |(|). Note, that on the 


contrary to the QCD /3-function, analytically cal¬ 
culated in Ref. [ |l3) at the N 3 LO level, the ex- 
pression for 7 /./ is known up to NNLO order. 
Moreover, its NNLO corrections were calculated 
in case of odd n < 13 only. In order to fix the 
numerical values of the NNLO corrections to 7 p ' 
(and thus the term q(n)) for even n inside the 
interval 3 < n < 13 we used the smooth interpo¬ 
lation procedure, proposed in Ref. [ 0 , and sup¬ 
plemented it by fine-tuning of definite NNLO co¬ 
efficients of 7 p 1 ' 1 [ |(| . The application of this pro- 

(n) 

cedure for estimating NNLO coefficients of 7 ^ 
with n even and Nf = 4 result in the numbers, 
which differ from the corresponding NNLO terms 
of non-singlet contributions to 7 ^ [ |l5[ |l^] in 
the 4th significant digit. The NNLO correction 
to 7 p' was estimated by us using extrapolation 
procedure, which has definite theoretical uncer¬ 
tainties. As to the applicability of the smooth in¬ 
terpolation procedure, we checked that it is repro¬ 
ducing the known even contributions to ( n ) 

with satisfactory precision [ [|. That is why we 
consider the results of its application, including 
the estimates of even terms of C^(n), as really 
reliable. 

Fixing by this way the N 3 LO coefficients 
C^ 3 \n) and estimating N 3 LO correction r{n ) to 
Eq. ( 8 ) by means of [1/1] Pade approximation 
technique, previously used in perturbative QCD 
e.g. in Ref. [ [ilj], we can use Eq. (1) for perform¬ 
ing approximate N 3 LO fits to X.F 3 data. At the 
next page the results, obtained in Ref. [ in the 
case of combining Eq. (1) with the IRR model of 
Eq. (3), are presented for N max = 6, first studied 
in Refs. [ f|, f|, and N max = 9 (see Table 1). 

Looking at Table 1 we arrive at the following 
conclusions: 

1 ) The NLO fits seem to support the IRR model 
of Ref. [ [ij by the foundation of the negative 
values of A 2 . 

These values are in agreement with the results 
of the previous fits of Refs. [ [|, P and with 
the ones, obtained in Ref. [ |17| using the NLO 
DGLAP analysis of the same set of CCFR’97 data 
and the parton distributions set (PDFs) of Ref. [ 
0 . The similar value A 2 = —0.104±0.005 GeV 2 
was also found in the NLO fits to the combined 
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Table 1 

The results to the fits of CCFR’97 data. A 2 [GeV 2 ] is the IRR model parameter. A^IL is measured in 
MeV. The cases of different Q g and N max are considered. 


order /N max 

O 

0 to 

II 

5 GeV 2 

20 GeV 2 

100 GeV 2 

NLO/6 

Ai^_ 

MS 

X 2 /nep 

A 

370±38 

80.2/86 

-0.121i0.052 

369T41 
80.4/86 
—0.121i0.053 

367i38 
79.9/86 
—0.120T0.052 

NLO/9 

A HP 

MS 

X 2 /nep 

A 

379i41 
78.6/86 
—0.125T0.053 

376T39 
79.5/86 
—0.125i0.053 

374i42 
79.0/86 
—0.124i0.053 

NNLO/6 

aM 

MS 

X 2 /nep 

A 

297i30 
77.9/86 
—0.007±0.051 

328T36 
76.8/86 
—0.017i0.051 

328i35 
79.5/86 
—0.015T0.053 

NNLO/9 

A W_ 

MS 

X 2 /nep 

a 2 

331i33 
73.1/86 
—0.013i0.051 

332T35 
75.7/86 
—0.015i0.051 

331i35 
76.9/86 
—0.016±0.051 

N 3 LO/6 

a^L 

X 2 /neP 

A 

305i29 

76.0/86 

0.036i0.051 

327T34 

76.2/86 

0.033i0.052 

326i34 

78.5/86 

0.029T0.052 

N 3 LO/9 

A SL 

MS 

X 2 /nep 

A 

333i34 

73.8/86 

0.038i0.052 

328T33 

75.9/86 

0.035i0.052 

328i38 

76.4/86 

0.034i0.052 


F 2 data [ [L9) using MRS(R2) PDFs [@. 

2) At the NNLO the values of A 2 are compara¬ 
ble wit zero within statistical error bars. The sim¬ 
ilar small value, namely A 2 = —0.0065 ± 0.0059, 
was obtained from the NNLO fits to F 2 data in 
Ref. [||. 

3) The inclusion of the N 3 LO corrections make 
A 2 positive. However, it has the statistical uncer¬ 
tainties twice as large as the central value. 

Thus we conclude that starting from the NNLO 
the IRR-model corrections to xF 3 can not be ex¬ 
tracted from CCFR’97 data with reasonable pre¬ 
cision and are shadowed by perturbative QCD ef¬ 
fects. 

4) The values of y 2 decrease from NLO up to 
NNLO and at the N 3 LO it almost coincide with 
the ones obtained at the NNLO. Moreover, y 2 de¬ 
creases with the increase of N max . This is the wel¬ 
come feature of including in the fits more detailed 
information on the perturbative theory contribu¬ 
tions both to the coefficient functions and the 


anomalous dimensions of xF 3 moments. 

5) For N max = 9 AM and A 2 are rather sta¬ 
ble to variation of Qq not only at the NLO but 
at the NNLO and N 3 LO as well. This property 
gives favor of our new results from Ref. [ |(ij| in 
comparison with the ones obtained in Ref. [ |i| 
for N max = 6 and Q\ = 20 GeV 2 using more 

/Q\ 

approximate model for 7 ( n ) and Pade approx¬ 
imation for Cp^(n). 

To transform A-^L into the values of a s (Mz) 

we first used the MS- scheme matching condition 
of Ref. [ pT) with the matching point chosen 
as m 2 < < 36to 2 following the proposal of 

Ref. [ ^2|. This gives us the possibility to esti¬ 
mate threshold uncertainties in a s (Mz)- Vary¬ 
ing the factorization and renormalization scales 
Fr = Ffac = A in the interval /xf^/4 < < 

4/x Ws we estimated the scale-dependent uncer¬ 
tainties. As the result we obtained the following 
values of a s {Mz) extracted from the fits to xF 3 
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CCFR’97 data with 1/Q 2 non-perturbative cor¬ 
rections modeled using the IRR approach [ [| : 

NLO a s (M z ) = 0.120 ± 0.002 (. stat) 
±0.005 ( syst) ± 0.002 (thresh)t.om6 (scale) 
NNLO a s (M z ) = 0.119 ± 0.002 (stat) (9) 
±0.005 (syst) ± 0.002 (thresh)^ 0^02 (scale) 
N 3 LO a s (M z ) = 0.119 ± 0.002 (stat) 
±0.005 (syst) ± 0.002 (thresh)^' qq 2 (scale) 


The systematic errors are fixed from separate con¬ 
sideration of these experimental uncertainties of 
the CCFR’97 collaboration. Notice, that the in¬ 
clusion of higher-order perturbative QCD correc¬ 
tions into Eq. (6) minimize essentially the scale- 
dependence uncertainties. They are in agreement 
with the similar estimates, obtained in Ref. [ |23j 
in the process of the fits to the definite model 
of xF 3 data. Our NNLO value of Eq. (10), ob¬ 
tained in Ref. [ |cj , within existing error-bars is in 
agreement with the results of NNLO Bernstein 
polynomial fits to the CCFR’97 data, namely 
a s (M z ) = 0.1153 ±0.0041 (exp) ± 0.0061 (theor) 

[| and a s (M z ) = 0.1196™! [|- 

In Table 2 we present the outcomes of our NLO 
and NNLO fits to the subset of CCFR’97 data, 
analyzed in Ref. [ ||. The twist-4 term was fixed 
with the help of the gradient model of Eq. (4). 


Table 2 

The results of the fits to the subset of CCFR’97 
data with HT defined by the gradient model with 
the coefficient B 2 [ GeV 2 ]. 


order 

(MeV) 

B' 2 ( HT) 

X 2 /nep 

NLO 

371±72 

-0.135±0.113 

75/74 

NNLO 

316±51 

-0.031±0.088 

64/74 


At the NLO the value of B 2 is in agreement with 
the value of the IRR model parameter A 2 . At the 
NNLO B 2 is comparable with zero. Thus we con¬ 
firm the existence of the effect of the shadowing 
of the dynamical l/Q 2 -corrections to xF 3 SF at 
the NNLO of perturbative QCD. This effect was 
first observed using the IRR model in Ref. [ [|] . 

If we use model-independent parameterization 
of the twist-4 contribution (see Eq. (5)), this ef¬ 
fect is becoming even more vivid. The results of 
extraction of h(x) in different orders of perturba¬ 
tive QCD and for different N max are presented at 


Fig.l, taken from Ref. [ ^J. 


0 2 = 20 GeV 2 
0 
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Figure 1. The x-shape of h(x) extracted from the 
fits of CCFR’97 for Q% = 20 GeV 2 . 

Several comments are in order. 

1) The x-shape of h(x), obtained at the LO 
and NLO, is in satisfactory agreement with the 
prediction of the IRR model of Ref. [ Q] . 

2) In all orders of perturbation theory the re¬ 
sults are rather stable to the variation of N max . 

3) The x-shape of h(x), obtained during the 
NNLO and approximate N 3 LO fits, demonstrate 
oscillation-type behavior with large error-bars. 
Thus we conclude, that starting from the NNLO 
the x-shape of h(x) is strongly correlated with 
higher-order perturbative QCD corrections to Eq. 
( 8 ). 

However, it is possible that more detailed un¬ 
derstanding of the NNLO behavior of h(x)/Q 2 
contribution to xF 3 will be obtained after NNLO 
analysis with taking into account systematic un¬ 
certainties of the data [ 241. 
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To conclude, we demonstrated that the in¬ 
clusion into the fits to CCFR’97 xF 3 data of 
the NNLO perturbative QCD effects is leading 
to effective shadowing of the dynamical 1 /Q 2 - 
corrections. However, it might be possible, that 
they will be detected in future even at the NNLO. 
This might happen in case more precise experi¬ 
mental data for the SFs of vN DIS will be ob¬ 
tained, say at the future neutrino factories (for 
detailed discussions see Ref. [ |25|). 

Note also, that while considering massive- 
dependent perturbative series for the Adler func¬ 
tion of e + e“-scattering [ |2(| and re-extracting 
the gluon condensate value from the charmonium 
sum rules [ |2lj with the calculated in Ref. [ ^8) 
three-loop massive corrections to their spectral 
function the effects of influence of higher-order 
perturbative QCD effects to the value of the gluon 
condensate, introduced in Ref. [ ||i|, were ob¬ 
served. Thus, the problems of correlations be¬ 
tween perturbative and non-perturbative terms, 
discussed in our talk, seem to be typical to other 
cases also, when it is necessary to analise the ex¬ 
pansions in both a s and inverse powers of Q. The 
similar point of view was expressed at this Sym¬ 
posium in the talk of Ref. [ [f0| . 
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